ACEEE Int. J. on Signal & Image Processing, Vol. 03, No. 01, Jan 2012 



IFKSA-ESPRIT - Estimating the Direction of Arrival 
under the Element Failures in a Uniform Linear 

Antenna Array 

Yerriswamy T. 1 and S. N. Jagadeesha 2 
'P.D. Institute of Technology, Visvesvaraya Technological University, Belgaum, Karnataka, India 

swamy_ty@yahoo.com 

2 J. N. N. College of Engineering, Visvesvaraya Technological University, Belgaum, Karnataka, India 

jagadeesha_2003 @yahoo.co.in 



Abstract — This paper presents the use of Inverse Free Krylov 
Subspace Algorithm (IFKSA) with Estimation of Signal 
Parameters via Rotational Invariance Technique (ESPRIT) 
for the Direction-of-Arrival (DOA) estimation under element 
failure in a Uniform Linear Antenna Array (ULA). Failure of 
a few elements results in sparse signal space. IFKSA algorithm 
is an iterative algorithm to find the dominant eigenvalues 
and is applied for decomposition of sparse signal space, into 
signal subspace and the noise subspace. The ESPRIT is later 
used to estimate the DO As. The performance of the algorithm 
is evaluated for various elements failure scenarios and noise 
levels, and the results are compared with ESPRIT and Cramer 
Rao Lower bound (CRLB). The results indicate a better 
performance of the IFKSA-ESPRIT based DOA estimation 
scheme under different antenna failure scenarios, and noise 
levels. 

Index Term s — Inverse Free Krylov Subspace, ESPRIT, 
Direction-of-Arrival. 

I. Introduction 

Antenna Array Signal Processing and estimation of 
Direction-of-Arrival of the signals impinging on the array of 
sensors or antennas, is a major functional requirement in 
radar, sonar and wireless radio communication systems. 
Generally, a Uniform Linear Antenna array (ULA) of a few 
tens to a few hundred elements (large antenna array) [1] are 
used for processing the signals impinging on the antenna 
array and estimate the DOAs. Among the various high 
resolution methods for DOA estimation, subspace based 
methods are most popular and powerful method. The 
popularity is due to its strong mathematical model to illustrate 
the underlying data model and it can with stand the 
perturbations in the data [2]. Subspace methods for DOA 
estimation searches for the steering vector associated with 
the directions of the signals of interest that are orthogonal to 
the noise subspace and are contained in the signal subspace. 
Once the signal subspace is extracted the DOAs are estimated. 
The decomposition is performed using the Eigen Value 
Decomposition (EVD) of the estimated received signal 
correlation matrix. The Multiple Signal Classification (MUSIC) 
[3] and Estimation of Signal Parameters via Rotational 
Invariance Technique (ESPRIT) [4] are the popular subspace 
based DOA estimation algorithms. 
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MUSIC algorithm is a spectral search algorithm and requires 
the knowledge of the array manifold for stringent array 
calibration requirement. This is normally an expensive and 
time consuming task. Furthermore, the spectral based methods 
require exhaustive search through the steering vector to find 
the location of the power spectral peaks and estimate the 
DOAs. ESPRIT overcomes these problems by exploiting the 
shift invariance property of the array. The algorithm reduces 
computational and the storage requirement. Unlike MUSIC, 
ESPRIT does not require the knowledge of the array manifold 
for stringent array calibration. There are number of variants 
and modification of ESPRIT algorithm [5] [6] [7] [8] [9] . ESPRIT 
algorithm is also extended for sparse linear antenna arrays or 
non-linear antenna arrays [10] [1 1]. For nonlinear arrays the 
aperture extension and disambiguation is achieved by 
configuring the array geometry as dual size spatial invariance 
array geometry [10] or by representing the array as Virtual 
ULA, and using the Expectation-Maximization algorithm [11]. 
The subspace algorithms are heavily dependent on the 
structure of the correlation matrix and are unsuitable to handle 
sensor failures. 

For handling sensor failure many modifications are 
proposed. Larson and Stocia [12], proposed a technique for 
estimating the correlation matrix of the incomplete data using 
the ML approach and shown improvement in MUSIC for 
handling sensor failure. However, it increases the complexity. 
A method for DOA estimator to handle sensor failure based 
on neural network is proposed by Vigneshawaran etal [13]. 
The technique can handle correlated signal sources, avoids 
the eigen decomposition. The drawback with these 
techniques is initialization of the network, and is performed 
by trial and error method. The authors in [15] proposed a 
DOA estimation technique by combining the EM algorithm 
with the MP method. The EM algorithm, expects the missing 
data and maximization is performance using the MP method. 
The method suffers for the drawback of increased complexity 
and requires good initialization. In [20], authors have 
proposed a Direct Data Domain DOA estimation algorithm, 
using the matrix completion procedure and Matrix Pencil (MP) 
method to first estimate the missing data observations and 
then estimating the DOA. However, the MP method is highly 
sensitive to the perturbations. Modification of ESPRIT DOA 
estimator to handle sensor failures is not present. 
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Decomposing the sparse correlation matrix is required to make 
the ESPRIT able to estimate the DOA for a faulty or incomplete 
data resulted from the failure of sensors. 



incident wsvetront 




Reference 
Element 

Figure 1 : Illustration of ULA for DOA estimation 

The very popular technique to solve the eigenvalue problem 
of a large sparse matrix is known as Krylov Subspace based 
techniques [15]. The reason for its popularity is due to its 
generality, simplicity and storage requirements. The iterative 
techniques of Krylov Subspace based techniques include 
Arnoldi approach, Lanczos algorithm, Jacobi-Davidson 
algorithm and Inverse Free Krylov Subspace Algorithms 
(IFKSA) proposed by Qin and Golub [16]. The study and 
analysis of the latter three techniques can be found in [15] 
[17]. The Krylov Subspace methods are applied in wide areas 
of application. For example, the method in used in control 
systems for estimating the parameters of a large sparse matrix, 
electric networks, segmentation in image processing and 
many others. IFKSA iteratively improves the approximate 
eigen pair, using either Lanczos or the Arnoldi iterations at 
each step through Rayleigh-Ritz projection procedure [17]. 
The algorithm is a very attractive due to the following 
reasons; first, the technique can be used to find any number 
of smallest eigenvalues (Largest can also be calculated), and 
second, the algorithm is less sensitive to perturbations. 

The main objective of this paper is to extended 
conventional ESPRIT algorithm to handle the element failures. 
The technique uses IFKSA method to decompose the 
estimated correlation matrix in to signal subspace which 
corresponds for larger eigenvalues and using the signal 
subspace the conventional ESPRIT is followed, resulting in 
IFKSA-ESPRIT, a robust DOA estimation algorithm to handle 
sensor failures. 

Before evaluating the IFKSA-ESPRIT algorithm for 
antenna failures, we first carried out a detailed evaluation of 
the performance of IFKSA-ESPRIT DOA estimator by 
comparing with the conventional ESPRT under different 
Signal-to-Noise Ratio (SNR) conditions and compared with 
Cramer Rao Lower Bound (CRLB) [18]. Next, the performance 
of IFKSA-ESPRIT scheme has been evaluated under different 
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array sensor failure scenarios and noise levels and the results 
are compared with CRLB. Finally, evaluating the performance 
of algorithm for large antennas, wide and narrow angles of 
arrival for various element failure scenarios and noise levels 
is performed. Results indicate that proposed technique has 
smaller errors for all the scenarios studied. 

The remainder of the paper is organized as follows. In the 
following section the signal model is discussed, followed by 
overview for Krylov subspace in section 3. In Section 4, we 
discussed the proposed IFKSA-ESPRIT algorithm followed 
by the simulation results in section 5. Finally, conclusions 
are discussed in section 6. 

II. Signal model 

The DOA estimation problem is to estimate the directions 
of plane wave incident on the antenna array. The problem 
can be looked as parameter estimation. We here mainly 
introduce the model of a DOA estimator. Consider an If - 
element uniformly spaced linear array. The array elements are 
equally spaced by a distance d, and a plane wave arrives at 
the array from a direction & off the array broadside. The angle 
is called the direction-of-arrival (DOA) or angle-of-arrival 
(AOA) of the received signal, and is measured clockwise 
from the broadside of the array. 

Let N narrowband signals all centered around a known 
frequency, impinging on the array with a DOA , . The received 
signal at the array is a superposition of all the impinging 
signal and noise. Therefore, the input data vector may be 
expressed as 

w 



: (0 =]T a <0 i > i <O+a(£} 



Where, 



e -]^-d(M--Osin9~\ 



(1) 



(2) 



a(f?i) is the steering vector of the antenna array. Here, T 
represents the complex conjugate transpose. 

X = A(8)s + W (3) 

Where , X is the array output matrix of size M X K, A is the 
complete steering matrix of size M X N function of the DOA 
vector 3, s is signal vector of size ± ty x K and W is the noise 
vector of size . Here, is the number of snapshots. Eq. (3) 
represents the most commonly used narrowband input data 
model. When an element fails there will be output from the 
failed element. The received sparse signal vector due to failure 
of elements is 

X' = A'(9)s + W (4) 

Where, If" is number of elements functioning. The correlation 
matrix of the sparse received signal is 

K 

1 V 

(5) 
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The objective of the papers it to estimate the DOA means to 
find out the value of the parameter , when a few of the 
elements fail to work. 

III. Overview of Krylov Subspace Method 

We use Krylov subspace method [15] [17] to decompose 
the received signal space in to signal subspace and noise 

Table I. Inverse free preconditioned krylov subspace - esprit algorithm 
(ifks-esprit) 

Input: R'. Xuniber c: SLsnsl Sources 

Stepl: Compute the Eigen decomposition of me array covariauce 
matrix using the Inverse Free preconditioned Krylov 
subspace (IFKS) algorithm 

IFHAlgorifonfJ: 

Step a: Set m > 1, * <- R'and an initial approximation 

u„ with |Kl| = i 

Step b: F or k = 0, 1, ■ ■ ■ ., tmti H Coni ergsnce 

C onstruct a basis V m = [v a , v 1 i>J for the 
subspace 

JC* = sptm(u^ (* - 1^)1%. (* 

Stepc: Form* m = V m T (*-pgv„ 
Step d: Compute the {i F j,',) eigenpairs of <P m . and select 
the desired ones. 

Pfc-l = Pk + K Hk-l = VmK 

End: 

Step 2: Estimate the D OA usins: ESPRIT 
ESPSITf): 

Step a: Fcrm the signal subspace E using the 

eigenvectors obtained from step 1 
Step b: Form the matrices E„ and E t 
Step c: Estimate the D OA solving [E„, E : ]. 
End: 

subspace spanned by the dominant eigenpairs and the 
smaller eigenpairs respectively. The decomposing is the first 
step in ESPRIT DOA estimator. A Krylov subspace of 
dimension n. generated by a vector u and the matrix <1> is, 

JC n (u, $) = span(a, $ 2 a, ■ ■ ■ , £ n_1 u) (6) 

Where, uis the initial eigenvector. The idea is to generate an 
orthonormal basis V = (v L v ; .. — .. v n }, of 3C ™ and seek an 
approximation solution to the original problem from the 
subspace. Since we seek approximate eigenvector 3, u E Jf , 
we can write 

u = Vy (7) 

The approximate eigenpair (Xu) is obtained by imposing 
Galerkin condition [17] 

vf (SVy- AVy) = 0, v< E X j = %2 — ,N (8) 
Or 

(®-M)u±K« (9) 

Therefore y and i. must satisfy 

B,j=Ay (10) 

Where B„ = V d *V 

The approximate eigenvalues i. resulting from the projection 
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process are all the eigenvalues of the matrix EL- The 
associated eigenvectors are the vectors Vy in which y is an 
eigenvector of EL associated with the i. . The procedure for 
numerically computing the Galerkin condition to eigenvectors 
and eigenvalues of is known as Rayleigh-Ritz procedure 
[15]. The approximation is updated convergence. There are 
many algorithms used to solve the above problem, namely 
Arnoldi algorithm, Lanczos algorithm and Jacobi-Davidson 
algorithm. Golub and Yi [8] proposed another algorithm known 
Inverse Free Krylov Subspace Algorithm. Inverse Free Krylov 
Subspace method starts with the initial approximation and 
aims to improve the approximation by minimizing the Rayleigh 
Quotient on that subspace. 

IV. LFKSA-ESPRIT Algorithm 

The algorithm for estimating the DOA from a faulty ULA 
is given in Table 1. It consists of two steps. Step 1 is 
decomposing the correlation matrix for signal subspace and 
the noise subspace using IFKSA technique. In the second 
step, the conventional ESPRIT algorithm in applied to 
estimate the DOAs. The algorithm starts with an initial 
approximation {p B ,u D ) and aims it improving through 
Rayleigh-Ritz projection [16] on a certain subspace, i.e., by 
minimizing the Rayleigh quotient on that subspace 



p(u) = 



u r $u 



(11) 



The Rayleigh quotients are important both for theoretical 
and practical purposes. The set of all possible Rayleigh 
quotients as uruns over C is called the field values of The 
gradient of Rayleigh quotient at u p is 



Vp(u ) = p D = V 



(u fl r «u a )u a r u a - n i r *n,g— u, r u, 



2((3>Uo) -p u ) 



(12) 
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The well known steepest descent method chooses a new 
approximate eigenvector u , E 5jjtm[u D , Opjby minimizing 
pCtiiX Clearly, this is a Rayleigh-Ritz projection method on 
the subspace 3f> = spcm[u 0l (# - p D )u } The 
new approximation u L is constructed by inexact inverse 
iteration [9] 

u, =($-p D )- 1 u D (13) 
then u- is indeed chosen from a Krylov subspace as generated 
by (<|> — pjj. Therefore, a new approximate eigenvector is 
found from the Krylov subspace of some fixed n 
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0C 71 = span^iiQ, ($ 



(14) 



(18) 



- P ) n 

by using the Rayleigh-Ritz projection method. The improved 
approximate eigenvector is obtained by iterating the above 
steps. The step 2 of the algorithms is the conventional ESPRIT 
algorithm. ESPRIT [4] is a computationally efficient of DO A 
estimation. In this technique two sub-arrays are formed with 
the identical displacement vector, that is, in same distance 
and same direction relative to first element. The array 
geometry should be such that the elements could be selected 
to have this property. The ESPRIT 




Fig. 1: RMSE vs SNR plot of CRLB, ESPRIT and IFKSA-ESPRIT 
algorithms for M = 1 1 and all elements are working 

algorithm uses the structure of the ULA steering vectors in a 
slightly different way. The observation here is that A has a 
so called shift structure. Define the sub-matrices, A - and A . 
by deleting the first and last columns from respectively, and 
note that and are related as 

A 1= A J1 (15) 

Where, £1 is a diagonal matrix having the "roots" e-?*f, 
j = Li,— If. Thus DO A estimation problem can be reduced 
to that of finding . 

Let and are the two matrices obtained by deleting the 
first column and last column of (5). Also let and are the two 
matrices with their columns denoting the eigenvectors 
corresponding to the largest eigenvalues of the two 
autocorrelation matrixes and respectively. As these two sets 
of eigenvectors span the same - dimensional signal space, it 
follows that these two matrixes and and are related by a 
unique nonsingular transformation matrix , that is 

E 1 = E D T (16) 

These two matrices E Ci and E, are related to the steering 
vector A c and A L by another unique nonsingular 
transformation matrix r, as the same signal subspace is 
spanned by these steering vectors. Thus 



A r 



v 1 = = A nr 

Substituting for V ft and V L , one obtains 
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Which states that the eigenvalues of 1* are equal to the 
diagonal elements of .f£ and that the columns of f are 
eigenvectors of . This is the main relationship in the 
development of ESPRIT. 

IV. Simulation Results 

In this section we examine the performance of the 
proposed algorithm for estimating the DOAs, when a few 
elements fail to work. We consider a ULA of M elements, with 
inter-element spacing of. The 300 snapshots of the signals, 
which are assumed as complex exponential sequences 

s (t) = e ma t = 1,2, ...,300 (19) 

is considered for the simulation. The Signal-to-Noise Ratio 
(SNR) is defined as 



SNR = 10log 10 ^ 



(20) 



Where, cr„ : and rcjjj are signal power and noise power respec- 
tively. The Root Mean Square Error (RMSE) is used to evalu- 
ate the performance of the algorithm. 




Fig. 2: RMSE vs SNR plot for estimation of DOA for various 
location of failure of two elements by IFKSA-ESPRIT algorithm 
for M =11 elements 

The RMSE is defined as 
J 



RMSE = 



(21) 
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Where, f(, ) is the mean value, & is the actual DOAs vector 
and pis the vector of estimated DOAs. We considered three 
simulation examples. 

First, we evaluate the proposed method for the case when 
all the elements are functioning. The results are compared 
with the conventional ESPRIT and Cramer-Rao Lower Bound 
(CRLB). The SNR vs RMSE plot is plotted for evaluation of 
the algorithms. Two signals of equal amplitudes are assumed 
to be impinging in the array of size 1 1 elements for the DOAs 
[5° -5°]. It can be observed form Fig. (1), the algorithm is 
showing improved performance when compared to ESPRIT 
at low SNR values and also able to achieve the CRLB [ 10] . 

In our second example we evaluate only the proposed 
method. The CRLB is calculated for the failed elements. We 
assume to elements are failed at some fixed positions for 
calculating the CRLB for various SNR values. The proposed 
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algorithm is evaluated for failure of two elements at various 
locations. For example, the locations we selected are (3 rcl and 
7 th ) elements at the ends of the array, (1 st and 2 nd ) first two 
elements, (10 th and 1 1 th ) last two elements, and (5 th and 6 th ) 
middle two elements. With the failure of two elements the 
algorithm is able to estimate the DOAs and achieve CRLB, 
shown in Fig. (2). Furthermore, when the (3 rd and 7 th ) ele- 
ments fail, at low SNR values the algorithm shows perfor- 
mance degradation. 

In our final example, we consider a large array of size 100 
elements. The number of signals impinging on the array is 
taken as 10 of equal amplitudes. The actual DOAs of the 
signals are assumed as [-30° ,-20° ,-10° ,-5° ,0° ,5° ,10° ,20° ,30° 
,40°]. The algorithm performance for this large array is 
evaluated for 6 cases; all the 100 elements are functioning, 4 
elements are failed at the end i.e, at locations 97, 98, 99 and 
100, 4 elements are failed at the beginning i.e, at locations 1 , 
2, 3 and 4, 10 elements at random locations, 20 elements at 
random locations and finally 40 elements at random locations. 
The SNR vs RMSE plot is plotted in Fig. (3). It is observed 
that except the last case, of failure of 40 elements at random 
locations, the proposed algorithm is able to estimate the 
DOAs. When 40 elements fail to work the performance gets 
worse. 
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IV. Conclusion 

In this paper, the performance of the IFKSA-ESPRIT DOA 
estimator is evaluated for various antenna array element 
failures in noisy environment. Results indicate that the 
algorithm performance better than the conventional ESPRIT 
at low SNR values, when all the elements are functioning. For 
failure case, IFKSA-ESPRIT is able to estimate the DOAs. 
However, when there is large number of element failures the 
algorithm is failed to estimate the DOA, which need further 
investigation. 
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